Energy spectrum of the low-lying gluon excitations in the Coulomb gauge 



O 

o 

Oh 

<: 



> 

oo 

o 

o 
\o 

O; 

Oh. 
I 

a.: 



X 



Adam P. Szczepaniak and Pawel Krupinski 

Physics Department and Nuclear Theory Center 

Indiana University, Bloomington, Indiana 4^74-05 

(Dated: February 2, 2008) 

We compute the energy spectrum of low-lying gluonic excitations in the presence of static quark- 
antiquark sources using Coulomb gauge and the quasi-particle representation. Within the valence 
sector of the Fock space we reproduce both, the overall normalization and the ordering of the spin- 
parity multiplets. We discus how the interactions induced by the non-abelian Coulomb kernel are 
central in to fine structure of the spectrum. 
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I. INTRODUCTION 

Lattice simulations of the energy of static quark-quark 
or quark-antiquark {QQ) systems have been studied for 
various configurations of the gluonic field [E S S 111- 
They effectively represent the spectrum of gluonic exci- 
tations in the presence of sources. The results of simu- 
lations indicate that the spectrum is quite complex with 
an unusual ordering of levels and unexpected behavior for 
large separations between the sources. In the following 
we will specialize to the spectrum in the presence of QQ 
sources and in the absence of light quark pairs. The wave 
function of the gluon field can be characterized accord- 
ing to its behavior under the symmetries of the QQ sys- 
tem. These include rotation around the QQ axis, hence 
after chosen to be in the z direction. The correspond- 
ing conserved quantum number, A represents projection 
of the total gluon angular momentum on the QQ axis 
and is conventionally denoted by A = S,n, A, • • • cor- 
responding to projections, A = 0, 1, 2, • • • , respectively. 
Other symmetries are the combined product of parity 
and charge conjugation, PC, with eigenvalues denoted 
by g and u for PC = -1-1 and —1, respectively and re- 
flection in the x — z plane, Y = ±1. The representation 
space of QCD eigenstates in the presence of static QQ 
sources is then given by states labeled as \R, h.y'^) where 
R is the distance between the sources. The corresponding 
energies will be denoted by V^pc{R). 

From lattice simulations the following characteristics 
of the the spectrum emerge. As expected, the ground 
state has E+ quantum numbers and as a function of R 
is well described by the " Coulomb-|- linear" , Cornell-type 
potential. The first excited state has one unit of gluon 
angular momentum along the QQ axis and negative PC. 
The state with A = 1 and positive PC has energy higher 
by 400 MeV at i? = 1 fm. This fact alone can be used to 
discriminate between various pictures describing the dy- 
namics of gluonic modes. For example, if the gluon field 
is thought of as a localized quasi-particle that interacts 
with the quark sources |3 in a way that satisfies Casimir 
scaling 6], then it would be expected that the S wave, 
spherically symmetric gluon wave function has lower en- 
ergy than the P-wave with one unit of orbital angular mo- 



mentum. Since gluons have internal spin-parity quantum 
numbers J^^ = 1^^ this would result in Hg {PC ~ 1) 
state having lower energy than the n„ {PC = — 1) con- 
tradicting lattice results. A similar reversed ordering is 
observed in higher excitations as well. The Ag state has 
lower energy than A„ which in the constituent picture 
means that the gluon state with two units of orbital an- 
gular momentum has energy lower than the state with 
one unit. The same is true for excited S states. The 
S~ state which in the constituent picture has one unit 
of orbital angular momentum has lower energy than the 
first excited S'"^ with vanishing orbital angular momen- 
tum. Higher in the spectrum the E~ corresponding to 
two units of orbital angular momentum has lower energy 
than E+ with one unit of orbital angular momentum, al- 
beit the splittings between these higher excitations are 
smaller than for the 11 states. These inconsistencies be- 
tween lattice results and the constituent picture have 
been noticed in Q- The bag model Q seems to be 
doing better in this respect. The boundary conditions 
on the gluon field, which set it free inside the bag cav- 
ity make the TE mode with pseudo- vector, J^'~'' = 1^ , 
quantum numbers to have lower energy than the TM 
mode with J^*^ = 1 . This in turn leads to the en- 
ergy of n„ state to be lower than for the Hg state p. 
Finally in Q it was observed that the non-relativistic 
fiux tube model [13 also predicts the state with A = 1 
and PC = -|-1 quantum numbers to be the first excited 
state. This is because in the fiux tube model the gluon 
degrees of freedom moving in a plane transverse to the 
QQ axis have negative parity but, unlike the vector po- 
tential have positive internal charge conjugation. Thus 
both the bag model and the flux tube model give the right 
ordering of the spectrum of low lying gluonic excitations, 
albeit for different reasons. The constituent gluon pic- 
ture of \^ was based on the mean fleld representation of 
the Coulomb gauge QCD, however it did not take into 
account the interactions emerging from the non-abelian 
Coulomb potential. In |ll| it was shown how such inter- 
actions affect the E+ potentials. Here we extend that 
analysis to gluonic excitations with other symmetries. 
The main finding is that the non-abelian Coulomb po- 
tential is responsible for reversing the naive ordering ex- 



pected from two body quark-gluon interactions. We also 
find that the overall scale of gluonic excitations in the 
mean-field quasi-particle picture is somewhat higher but 
consistent with the lattice results thus making the con- 
stituent gluon model a viable representation of the low 
energy gluon dynamics. This representation also explains 
the degeneracies in the spectrum seen in the lattice data 
at small QQ separations. Such degeneracies are expected 
as the system becomes more spherically symmetric. At 
large QQ separations the gluonic wave function is ex- 
pected to be qualitatively different from that of a sin- 
gle or a few quasi-particle state [13, [131 • As separation 
between color sources increases the mean field Coulomb 
interaction is expected to rise more rapidly than the true 
energy of the system fM, 15, 16|. Thus as the separation 
increases states with a large number of gluons separated, 
on average, by a small fraction of the QQ distance R are 
expected to have lower energy than states with a small 
number of gluons separated by a distance of the order of 
the QQ separation. At what QQ distances the transi- 
tion from the constituent, few gluon picture to the flux 
tube or string like picture takes place is however still an 
open question. Analysis of lattice results for the splitting 
between gluonic levels at large R does not conclusively 
favor the string like picture even for QQ separations as 
large as a few fermi [ij, [l^ . In [ij we have shown that 
if the mean field Coulomb energy also rises linearly with 
the R, as indicated by lattice computations [lj| then it is 
difficult to generate the constituent string and the state 
is dominated by valence constituent gluons. This is be- 
cause the Coulomb energy of a state of Ug gluons sep- 
arated from each other by a distance 0{R/n) does not 
depend on the number of gluons [13, thus an increase 
in the average number of gluons can only originate from 
mixing. However as R increases it turns out that mixing 
between states with a different number of gluons is much 
smaller than the diagonal Coulomb energy [111. Lattice 
results of J^] show that the Coulomb energy rises lin- 
early with the QQ separation with a string tension which 
is roughly three times bigger than that of the true energy, 
however, since the lattice QQ state used to compute the 
Coulomb energy is not the same at the mean field state 
which defines the quasi particle gluon basis it is possi- 
ble that the Coulomb energy of the latter grows with R 
faster than linearly, i.e as R'^ with 7 > 1. In this case 
the average number of gluons at fixed R would roughly 
correspond to the minimum with respect to Ug of 



En,{R) ^Ugmg+Jlg ( — 



(1) 



leading to Ug ^ R and energy rising linearly with R. 
Here nig represents the average kinetic energy of a quasi- 
giuon. It is thus possible that there is a smooth connec- 
tion between the quasi-particle and flux tube represen- 
tation, but we leave a more quantitative description for 
the future. In the following we explore the flue struc- 
ture of the Coulomb interaction and its role in ordering 
the spectrum of low-lying gluonic excitations. We will 



therefore restrict the quasi-gluon Fock space to a the sin- 
gle particle sector and fit the expectation value of the 
mean field Coulomb interaction to the ground state QQ 
energy [ll[2i [23. 

The paper is organized as follows. In the next sec- 
tion we summarize the basics of the mean field Coulomb 
gauge approach and describe the relevant interactions. 
In Section HI we present out results for the spectrum. A 
summary and outlook are given in Section IV. 



II. COULOMB GAUGE QCD HAMILTONIAN 
IN THE QUASI-PARTICLE REPRESENTATION 

The derivation of QCD in the Coulomb gauge is given 
in i23|. The canonical approach is to start from v4°'"(x) — 
the Weyl gauge and using the residual gauge free- 
dom perform a coordinate transformation to coordinates 
A°(x) constrained to satisfy V • A'^(x) =0 and N^ — 1 
phases (^"(x), a — 1,- ■ ■ , N^ — 1. In the new coordinates 
Gauss's law can be used to eliminate the dependence on 
the gauge phases and in the Shrodinger representation 
the QCD spectrum is formally obtained by solving 



ij[A«(x),n''(x)]^„[A«] = £;„*„[A^, 



(2) 



with the canonical momenta, 11" — —id/dA°'{x.), sat- 
isfying [Hf(x),A5(y)] = ^tSabS'^{V)dH^-y) where 

(5^(V) = (5*^— V^V-'/V^. The coordinate transformation 
from the Weyl gauge to the Coulomb gauge is nonlinear 
and leads to the Faddeev-Popov (FP) determinant in the 
space of the field configurations, J — Det{l — A), where, 
(1 - A)(x, a; y, b) = SabS^{^ - y) - (5/47r)/,,bV,(l/|x - 
y|)A'^(y). More discussion of the topological properties 
of the fundamental domain of the gauge variables can be 
found in [23 . The role of FP determinant has been in- 
vestigated in [ill I24L [23 where it was found that it can 
be effectively absorbed into the parametrization of the 
vacuum wave functional. In the following we will thus 
set J^ — I and use the ansatz for the variational (unnor- 
malized) vacuum wave functional of the form. 



(A|0) = *o[A] = cxp (^- I ^A'^{kMk)A^i-k)^ , 

(3) 
where A°(k) = Jdxexp(— zk • x)A"(x) is the Fourier 
transform for the coordinate space. 

The variational parameter a;(/c) is obtained by solv- 
ing the Dyson (gap) equation arising from minimizing 
the vacuum expectation value (vev) of the Hamiltonian 
d{0\H\0)/duj{k) = 0. The solution is well approximated 
by Lu{k) — nig for k = |k| < nrig and aj(fc) = k for k > nig 
with nig — 600 MeV |2l| . In computing the vev of the 



Hamiltonian the Coulomb energy 



He = iydxdyp^(x)AqA](x,a;y,6)p^(y), 



(4) 



contributes the energy of interaction between color 
charges Pg = —fabJ^''{^) • A'^(x). In Coulomb kernel, 
K[A\ appearing in He the self-interactions between the 
transverse gluons 

7^[A](x,a;y,6) = —j dz ^— ^ , (5) 

are evaluated to leading order in Nc- This enables to 
express the vev of K in terms of a set of two coupled 
Dyson equations whose solution can be written as 



\K[A]{^,a;y,b)\0)/{0\0)=-6abVc{\^-y) 



(6) 
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The expectation value of the kernel Vc is renormalized 
in such a way that the expectation value of the Coulomb 
interaction between QQ sources in the mean field gluon 
vacuum reproduces the lowest energy of the QQ state Ji 

The complete spectrum of gluon state can be obtained 
by successive application of quasi-particle, gluon creation 
operators, Q;°(k, A), defined with respect to the mean 
field vacuum, through the BCS transformation 



FIG. 1: Matrix elements, (i?,k', A'|H|_R,k, A). Diagrams a) 
and b) represent gluon and quark self energies, respectively. 
Diagrams c) and d) represent the Coulomb interaction, Vc 
between the gluon and one of the quarks and between the 
two quarks, respectively. In the bottom row, diagrams e) and 
f ) describe matrix elements of the interaction term resulting 
from expansion of the Coulomb kernel K[A] in up to one 
power in gluon field. 
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(27r)3 ^2^ 
-a''^t(-k,A)e(-k,A)] e 



ikx 



n"(x) 



dk 



'^^^^ fa"(k,A)e(k,A) 



(27r)3 V 2 
-a°'^(-k,A)e(-k,A)]e*'^'' 



(7) 



Here e represent helicity vectors with A = ±1. When 
describing the spectrum of gluons in the presence of QQ 
sources, we will truncate the quasi-particle gluon Fock 
space to contain a single quasi-gluon, i.e. the Hamilto- 
nian will be diagonalized in the basis of states spanned 
by 



|i?,k,A) 



1 



VN^f 



)t.,a^'"(k,A)r'^gl„, 



|o) 
(o|o)- 



(8) 



III. SPECTRUM OF LOW-LYING GLUONIC 
EXCITATIONS 



From the single quasi-particle state given in Eq. (jHl 
the states with good Ap^ quantum numbers can be con- 
structed 



for A 7^ and 



|i?,iV,0^^) = 



dk 

^2^ 



E i^'Nik)xix\R,K^) 



jg:?,M,A 



4^ ^Om W, 



(10) 



for A = 0. Here jg represents the total angular momen- 
tum of the quasi-gluon. It is a good quantum number 
only in the limit i? — > 0, while for finite R states with 
different values of jg can mix, although in our numeri- 
cal computations we have found that of a single jg state 
dominates the energy eigenstates for all values of R con- 
sidered. It is only the projection of the total angular 
momentum on the QQ axis. A, that is always conserved. 
The wave function xtx represents the two possibilities 
for the spin-oribt coupling of given parity. It is given by 



5^x/V2 for ^ = 1 Ug 



Lg±r 



and \Sf_i\/V2 for ^ = — 1 



{jg — Lg), corresponding to TM (natural parity) and TE 
(unnatural parity) gluons, respectively. The parity under 
reflection in the x — z plane, Y = ±1, for A 7^ is deter- 
mined by r]Y — ±1- The radial wave functions, ip''i^{k) 
are labeled by an excitation number N and jg and are 
solutions of the Coulomb gauge Hamiltonian projected 
onto the single quasi-gluon basis 



PHP\R,N,Alc) 



Vr 



N,At 



.iR)\R,N,Alc) 



(11) 



R,N,Alc)-^J^ E i^'Nik)xU\R,KX) 



2jg + 1 



Svr 



(27r) 



(9) 



Here P projects on the \qqg) states and Vj^j^pc{R) are 
the energies that will be compared with the la1;tice spec- 
trum. The matrix elements of PHP are shown in Fig. ^ 
and given explicitly in the Appendix. 

In terms of ^ and 77, the PC and Y quantum numbers 
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TABLE I: Spin-orbital wave functions of the lowest single 
quasi-gluon states. 



of the gluonic field are given by 



pc = a-^) 



Js + l 



(-1)^«, Y 



Cr/y(-l)^ for A 7^0 
C for A = 

(12) 
For A 7^ the two Y — ±1 states are degenerate. For 
small QQ separations the pattern of the spectrum mea- 
sured on he lattice can easily be understood since jg be- 
comes a good quantum number. In this case the gluon 
configurations of the eight lowest excitations are given 
in Tabled By setting i? = in the Hamiltonian matrix 
elements which involve gluon-quark interactions (see Ap- 
pendix) one finds the following. The interactions become 
A independent. This is expected since in the limit i? — > 
the system does not have a preferred direction while A se- 
lects one. The angular momentum barrier pushes states 
with higher orbital angular momentum, Lg up in en- 
ergy. The quark-gluon and antiquark-gluon interactions 
are attractive, while the quark-antiquark interaction in 
the color octet channel is repulsive. Thus the Coulomb 
tail of the quark-antiquark potential will eventually lead 
to rising energies as i? — > 0. Outside of this short dis- 
tance Coulomb region, however for given jg and ^ we 
expect states with different A to be degenerate. Further, 
the (^ — +1 multiplets (which contain Lg = jg — 1) are ex- 
pected to be lower in energy than the ^ = — 1 multiplets 
which have Lg = jg . This is indeed seen in the numerical 
results shown in Figs. I2I3I4I Compared to the lattice re- 
sults, however we see that while the degeneracy between 
different A states within ^-multiplets indeed occurs the 
order of the ^ = — 1 and ^ = +1 multiplets is reversed. 

To obtain the correct ordering it is necessary to include 
the irreducible, three-body interaction between quark, 
antiquark and gluon, shown in Fig. ^. This interac- 
tion is attractive and, as can be easily verified using the 
analytical expression in Eq. H2U|) , as i? ^ it pushes up 
in energy the ^ = — 1 {Lg ~ jg) states and has no effect 
on the £, = +1 states. As seen in Figs. 1516171 that show 
the results of diagonalization of the full Hamiltonian, this 
additional interaction energy is sufficient to change the 
order of the 11 and excited E states, and almost does the 
job for the higher energy, A states. 

When studying the large R behavior we first note that 
the long-range singularity of the linear potential is can- 
celed between the three self-energies and the two-body 
interactions. The three body interaction between the 
quark, antiquark and gluon is IR finite, which also im- 
plies that in the limit J? ^ oo it decreases with R. This 
is due to the gradient coupling of the transverse gluon to 
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FIG. 2: Spectrum of the A = 1, 11 states compared to the 
ground state QQ potential without the three body quark- 
antiquark gluon interaction corresponding to the diagram e 
in Fig.0 Lattice results are taken from y]. The energy scale 
is ro = (400 MeV)"^ 
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FIG. 3: Same as in Fig. |5|for the A = 0, E energies. 




FIG. 4: Same as in Fig.|21for the A = 2, A energies. 
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FIG. 5: Spectrum if the A = 1, 11 states compared to the 
ground state QQ potential using the complete Hamiltonian, 
which includes the three body quark-antiquark gluon interac- 
tions. 



FIG. 8: Splittings between the excited energies of quasi- 
gluons and the ground state QQ energy. The values A'^ ex- 
pected from the string model are, N = 1 for n„, N = 2 for 
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FIG. 6; Same as in Fig. |^for the A = 0, E energies. 
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FIG. 7: Same as in Fig. Ofor the A = 2, A energies. 



the Coulomb line. Furthermore, connecting one gluon to 
the Coulomb line effectively shortens the Coulomb po- 
tential making it less confining 11] . All this implies 
that splittings between excited states with Ug + 1 and 
Ug quasi-gluons at large QQ separations do not increase 
with R . They are expected to be roughly constant cor- 
responding to the average kinetic energy of a gluon in the 
color singlet state. Eg ^ uj{k ~ 1/^) = rUg ~ GOOMeV. 
On the other hand, to leading order in 1/R, separations 
of string excitations are expected to be proportional to 
AE = En+1 - En ^ Nn/R 0. In Fig. |H1 we plot 
{AE)/{Ntt/R) — 1 with AE representing the energy dif- 
ference between our Ug = 1 excited quasi-gluon energies 
and the ground state QQ energy. The corresponding val- 
ues of N were chosen as in [^ • The agreement with the 
lattice results, shown in Fig. 1 of [3| is very good. In the 
string model {AE)/{Ntt/R) — 1 is expected to approach 
zero at large separations, while the lattice and our results 
seem to indicate a positive slope. As discussed above this 
slope can be interpreted in terms of quasi-particle exci- 
tations. 



IV. SUMMARY AND OUTLOOK 

We computed the lowest excitation energies of the glu- 
onic field in the presence of static QQ sources. The 
gluons were described using single, quasi-particle states 
built above a mean-field variational ansatz for the ground 
state. The variational approach is expected to produce 
energy levels which are systematically higher than the 
true energies. Indeed the computed energies are above 
the lattice data by about one unit of Tq^ , i.e. 400 MeV. 
The nontrivial level ordering is however correctly repro- 
duced once the full effect of the Coulomb kernel is taken 
into account. In particular the irreducible three-body 



force is the one which is responsible of reversing the naive 
ordering expected in any constituent gluon model. The 
quasi-gluon gluon description also reproduces the degen- 
eracies in the energy levels for small separations, as well 
as the trends in the level splittings for intermediate sep- 
arations of the order of a few fm. At larger separations, 
multi-gluon effects are expected and those can be incor- 
porated into this picture by renormalizing the bare QQ 
potential. A quantitative analysis is currently being per- 
formed. 
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VI. APPENDIX 



The expansion coefBcients are computed by diagonaliz- 

mg Xne [JT^rnaxJ g^max ) -^ y^maxj g^max J niairiX, ^m'j';m,.jg^ 

obtained by evaluating the diagrams in Fig. Q] 



-^3 — H^a + H:n, + ■ ■ ■ + Hi 



3/: 



(15) 



evaluated in the basis of functions (j^m.j ■ In numerical 
computations for each jg, we used a momentum grid as 
the basis functions. The numerical results presented were 
for a single jg determined from Eq. (|12|l after verifying 
that increasing jg changes the computed spectrum by at 
most a few percent. For arbitrary Ap^^ the Hamiltonian 
matrix elements are given by 



H^a — 



^j'-Jc, f dkk 



.^—^^l^,^^^{k)Eg{k)<p,n,jJk)- (16) 



For each jg, the wave functions ipNj (k) are expanded 
in a complete orthonormal basis of functions (/>m j (k) 



The single gluon energy, Eg{k) is given in |llLl21| . 



i'Nj.ik) = Y^ a^j,'/>mjg(fc), 
m— 1 

with normalization, 



(13) 



Hsb = —CFVc{0)Sm' ,7nSj'^Jg 



= -AirCf 



dkk^ 

(2^ 



ycik)Sm',mSj' jg, (17) 



dkk^ 

(2^ 



■<P*m',j' (fc)'^mjg (k) = S,n',,nSj' Jg • (14) 



with Vc{R) fitted to the ground state QQ potential [I^, 

El. 






dq f dk 



2 ^ J (27r)3y (27r)3^™^J8 



</'m'j'(9)0m,ig(fc) / dx 



(2j'+l)(2jg + l) , . 



T/c(x-|) + Vfc(x+|) 



167r 



Nc V- f Jq_ f dk 

2 ..^„,J (2^)'i (27r)3--'.^^ 



C',,'(9)^c(k-q) 



g-jf.(k-q) _^gif.(k-q) 



^{2j'g + l){2jg + l) 

1677 



„ix-(k-q) 



^A,.'(q)^AA' (k)xL'XiA'^i;(q)^iA(k) + VYv'yi^ ^ -A) 



Xk) 



<.'(q)^AA' (k)xL'xiA'^i;(q)^iA(k) + vyv'y{^ - -a) 



ujiq) ^ V u;{k) 



/^(fc) , HA 

u;{q) ^ V u;{k) 



(18) 



and r]Y and ^ related to jg and Ap^ through Eq. H12|l . 



H: 



3d 



2N( 



-Vc{R)Srn',niSj' Jg 



1 f dkk^ 
= -^TT^TT^ I T:r-7^Vcik)joiRk)S,n' ,mSj, j^. 



2Nc J (27r)3 



(19) 



_ >r^ f dk dp dq fe'.j^ (p) (/)m.jg (fc) 



X / djcdydz 



"R R 

if (x - - , z + y - X, y + -) + (R ^ -R) 



gjxkgizqg-jyp 



(2j-+l)(2j, + l) 



Stt 



.«' ni 



<.'(p)^i';(p)xlv^i,o(q)^i%'(k)i^ia(k)xi'A^;;o(q) + '7y^y(A - -A) 



-E 



dk dp dq 'l>*m'.j',iP)^m.j,{k) 
(27r)3 (27r)3 (27r)3 ^2^^ ,/2^- 



K{k + q, q, p + q) [e'T-(k+P+2q) + (R ^ _R) 






^L'(p)^i:;(p)xi'.^i,o(q)^in*(k)^ia(k)xi,,<aq) +r7yry;.(A ^ -A) 



(20) 



r 



where the sum is over /z, i', A, A', cr, ct' and the kernel is and the kernel K(k,q,p) is given in |llLl2ll |. 
given by 



^^"'"'y) = /(^(^(^^(^''^'^) 



gixkgiypgizq 



(21) Finally, 
I 



X^ [ '^^ ^P ^"^ 

^■'^ " ^i 72^(2^72^ 



dk dp dq '/'m'j-(p)</)„,j,(fc) 



(2^)3 (2^)3 (27r)3 ^2^ TMfc) 
X / dxdydz 

y'(2j; + l)(2j, + l) 



if(x-5:,z + y-x,y-5:) + (R^_R) 



gixkgizqg-iyp 



Stt 



E 



<.'(p)^i:;(p)x^'.^i.o(q)^i:A'(k)i?i,A(k)xiA<'o(q) + ^i-^ylA 

if (k + q, q, p + q) [e^t '('^-P) + (R ^ -R) 



dk dp dq C'jv(p)</).^jjfc) 



(2^)3(2^)3(271)3 ^2^ TMfc)" 



y^(2j; + l)(2j, + l) 



Stt 



^AV'(p)^M;(p)xS'.^i.o(q)^i:A'(k)i?i.A(k)xi-A^i;o(q) + ^y^y(A - -A) 



« ni.*!- 



-A) 



(22) 
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